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A study of 0(l/mg) corrections for with lattice NRQCD 
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We investigate higher order effects in the nonrelativistic expansion of lattice QCD on the heavy-light meson 
decay constants and some other quantities in order to understand the truncation error of NRQCD. While our 
numerical results have large O(a) and 0(a s ) errors due to the use of Wilson light quark action and the tree-level 
matching, we find that the truncation error of higher order relativistic corrections are adequately small around 
the mass of the b quark. We also present a perturbatively matched results through 1-loop level without operator 
mixing effects. 



1. Introduction 

It is indispensable for the verification of the 
Standard Model and the new physics search that 
weak matrix elements between hadrons involving 
6-quark are calculated with high accuracy. Here, 
we focus on the heavy-light meson decay con- 
stants Q|, particularly Jb, whose accurate value is 
needed for the determination of |V^Vtd| together 
with B B parameter. The lattice NRQCD jg] en- 
able us to do the direct simulations at the b quark 
and has yielded remarkable progress on the heavy 
quark physics. So far the studies including 1/rriQ 
corrections have made it clear that there exists a 
large 1/toq correction to the value in the static 
limit Q. In this talk we report on the systematic 
study of the 1 /rriQ corrections. A preliminary re- 
port of an investigation of 0(1/toq) corrections 
similar to our work has been reported in Ref. 
The action and operators should be matched to 
the full theory since NRQCD is an effective the- 
ory. Recently we have completed the self-energy 
renormalization constants and the multiplicative 
part in the renormalization of axial vector cur- 
rent^]. In this article we present the results of 
/b and some other quantities with and without 
1-loop matching. 



2. Simulations 

Our numerical simulation is carried out with 
120 quenched configurations on a 16 3 x 32 lattice 
at /3 = 5.8. For the tadpole factor we employ 
uq = (Ppiaq) 1 ^ 4 with Ppiaq the average plaquette, 
which takes the value uq — 0.867994(13) mea- 
sured during our configuration generation. 

For the light quark we use the Wilson quark ac- 
tion with k=0.1570, 0.1585 and 0.1600. The chi- 
ral limit is reached at k c = 0.16346(7) and the in- 
verse lattice spacing determined from the rho me- 
son mass equals a -1 = 1.714(63) GeV. The hop- 
ping parameter k s corresponding to the strange 
quark is determined from m^/m p and rnx/m p , 
which yields k s = 0.15922(39) and 0.16016(23), 

respectively. We use the factor Jl — j^- as the 

field normalization for light quarkjfj]. 

For the heavy quark, we use the lattice NRQCD 
with the following evolution equation: 
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Where Hq is defined as follows: 
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and in order to realize two accuracies we define 
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two (5iJ's as follows: 
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The symbols A and A^ 2 ) denote the symmet- 
ric lattice differentiation in spatial directions and 
Laplacian, respectively, and A^ 4 ^ = J^ii^^) 2 ■ 
The field strengths B and E are generated from 
the standard clover-leaf operator. 

The original 4-component heavy quark spinor 
h is related to two 2-component spinors Q and 
X through Foldy-Wouthuysen-Tani(FWT) trans- 
formation, 



h(x) = R 



Q(x) 
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where R is an inverse FWT transformation ma- 
trix which has 4x4 spin and 3x3 color indices. 
We prepare two R's as before. After discretiza- 
tion, R at the tree level is written as follows: 
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where S- 1 — diag.-jV^ ', a^}. We apply the tadpole 
improvement to all link variables in the evolu- 
tion equation and R by rescaling the link variables 
as t/ M -> U^/uq. 

We performed two simulations using following 
two sets: 

set I = {6Hi,Ri} and set I = {SH 2 ,R. 2 }. (8) 

set I and I are fully consistent through 0(l/mg) 
and 0(1/ 17Iq), respectively. Furthermore 5H 2 has 
the leading relativistic correction to the disper- 
sion relation, which is an 0(l/mg) term, and the 
terms improving the discretization errors appear- 
ing in Hq and time evolution are also included. 
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Figure 1. 1/Mp dependence of fpMp 1 ^ 2 with set 
I(open symbols) and set I(solid symbols). Cir- 
cles are the results with tree-level matching and 
squares and triangles are with 1-loop matching on 
the scale of q*=Tr/a and 1/a, respectively. 

For the heavy quark mass and the stabilizing 
parameter, we use (arriQ, n)=(10.0,l), (5.0,2), 
(2.6,2), (2.1,3), (1.5,3), (1.2,3) and (0.9,4), which 
cover a mass range between 4m(, and m c . 

Our strategy is to study systematically the ef- 
fects of an inclusion of 0(l/mg) terms by com- 
paring two results from these simulations. All of 
our errors are estimated by a single elimination 
jack-knife procedure. 

3. Results 

3.1. decay constants 

Fig|l] shows the 1/Mp dependence of fp\JMp. 
Where Mp is a pseudoscalar heavy-light meson 
mass. Two vertical lines represent the mass re- 
gions corresponding to the B and D meson. Cir- 
cles (open for set I and filled for set E) are the re- 
sults with tree-level matching. From more quan- 
titative study, we can find that within tree-level 
analysis the relative magnitude of (l/m^) correc- 
tion is about 3%, conservatively 6%, around the 
B meson region |l. 

3.2. other quantities 

We also performed a similar investigation on 
Mp s — Mp and fpj fp, where P s denotes a pseu- 
doscalar heavy-strange meson. Those results are 
shown in FigJ| and ||. One can see from these 
figures that there is no significant difference be- 
tween the results from the two sets over almost 
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Figure 2. 1/Mp dependence of Mp B —Mp with set 
I(open symbols) and set II(solid symbols). Dia- 
monds are with k s from (f> meson and circles from 
K meson. Two asterisks indicate the experimen- 
tal values. 

all mass region. 

3.3. perturbative matching 

Recently we have completed the self-energy 
renormalization constants and the multiplicative 
part in the renormalization of axial vector current 
with set Since the scale q* where the strong 
coupling constant should be estimated is still un- 
determined, we estimate the matching factor at 
q* = 1/a and ir/a and show the results with the 
factor in Fig.|[ 

4. Summary 

We presented the effects of 0(1 /m^) correc- 
tion to the heavy-light meson decay constant with 
tree-level matched lattice NRQCD and Wilson 
quark action in a quenched approximation. While 
the 0(1/ mo) correction to the decay constant in 
the static limit is significant || , we find in our sys- 
tematic study that the 0(l/mg) correction is suf- 
ficiently small for B meson, so that there will be 
no need for incorporating 0(l/mg) corrections 
unless an accuracy of better than 5% is sought 
for. Our examination of other physical quanti- 
ties in the same respect also provides encouraging 
support to this statement. We have thus shown 
using our highly improved lattice NRQCD that 
the relativistic error, which has been one of the 
largest uncertainties in lattice calculations of the 
B meson decay constant, is well under control. 
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Figure 3. 1/Mp dependence of fpj fp with set 
I(open symbols) and set II(solid symbols). Dia- 
monds are with k s from (f> meson and circles from 
K meson. 

Our results still have several sources of large sys- 
tematic errors, that is, 0(a), 0(a s /M), 0(aa s ) 
and a quenching errors. When improvements of 
these errors are all in place, we expect to achieve 
a lattice NRQCD determination of fp with the 
accuracy of less than 10%. 

Numerical calculations have been done on 
Paragon XP/S at INSAM (Institute for Numeri- 
cal Simulations and Applied Mathematics) in Hi- 
roshima University. 
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